Spindown of Isolated Neutron Stars: 
Gravitational Waves or Magnetic Braking? 

Jan E. Staff 

Department of Physics and Astronomy, Louisiana State University, 202 Nicholson Hall, Tower 

Dr., Baton Rouge, LA 70803-4001, USA 

^N . Praslianth Jaikumar, Vincent Chan 

^vq ■ Department of Physics & Astronomy, California State University Long Beach, 

(^ . 1250 Bellflower Blvd., Long Beach, California 90840, USA 

^: 

l^ ■ Rachid Ouyed 

CN ■ Department of Physics and Astronomy, University of Calgary, 2500 University Drive NW, 

Calgary, Alberta, T2N 1N4 Canada 

• ■ ABSTRACT 

Of 

O ■ We study tfie spindown of isolated neutron stars from initially rapid rotation rates, 

-^ ■ driven by two factors: (i) gravitational wave emission due to r-modes and (ii) magnetic 

C^ ■ braking. In the context of isolated neutron stars, we present the first study including 
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self-consistently the magnetic damping of r-modes in the spin evolution. We track the 
spin evolution employing the RNS code, which accounts for the rotating structure of 



^^ . neutron stars for various equations of state. We find that, despite the strong damping 

o 



due to the magnetic field, r-modes alter the braking rate from pure magnetic braking 
for B < 10^'^G. For realistic values of the saturation amplitude asat, the r-mode can 

[■"^ ■ also decrease the time to reach the threshold central density for quark deconfinement. 

Within a phenomenological model, we assess the gravitational waveform that would 
result from r-mode driven spindown of a magnetized neutron star. To contrast with the 
persistent signal during the spindown phase, we also present a preliminary estimate of 

k>( , the transient gravitational wave signal from an explosive quark-hadron phase transition, 

\—i \ which can be a signal for the deconfinement of quarks inside neutron stars. 



1. Introduction 



Neutron stars are highly compact stars of typical radius i? ~ 12 km and mass M ~ I.SMq 
made mostly of degenerate neutron-rich matter at densities up to several times nuclear matter 
saturation density po = 2.5 x 10^^ g/cc. By tracking their long-term thern ial and rotationa l 
evolution, we can learn about the nature of matter under the crust. For example, iPage et al.l (|201ll ) 
have proposed that the thermal history of the neutron star in Cas A may be indicati ng the recent 
onset of neutron superfluidity deep in its interior. Recently, iNegreiros et al.l (1201 ll ) have shown 
how rotational evolution is linked to a reorganization of particle composition in the stellar interior, 
leading to switching on of exotic neutrino emission processes. A neutron star is a complex system 
with intertwined physical properties that can change on relatively short astrophysical timescales. 

In this paper, we will be concerned with the rotational evolution of a highly magnetized isolated 
neutron star. The electromagnetic emissions of a neutron star derive from its rotational kinetic 
energy, and its spindown is usually measured in terms of a braking index, n, which is dependent on 
the magnetic field configuration (n=3 for a dipolar fiel d). Neutron stars can also spindown through 
gravitational wave emissions associated to the r-mode (| AnderssonI 1 1 998l ) . In fact, the observational 
interest in r-modes comes from the fact that no neutron stars have been found to spin at rates near 
the maximum allowed frequency (the "break-up frequency"), r-modes offer a possible explanation 
of this fact: in rotating neutron stars, these modes lose energy through gravitational waves, which 
carry away angular momentum from the star and act as braking radiation. As the star spins down, 
its central density increases. This increase could be sufficient to make baryonic matter undergo 
phase transitions to more exotic phases of stroi igly interacting matter, such as quark niatter , 
with possible imp lications for gamma-ray bursts (JBerezhiani et al.l l2003l : iHaensel &: ZdunikI 120031 : 
Drago et al.l 120081 ) and the formation of quark stars, if t heoretical conjectures about the absolut e 
stability of strange quark matter are realized in nature (jltohl Il97d : iBodmeiJ Il97ll : IWittenI 11984 ) . 
While it is almost certain that this quark phase , if it e xists inside neutron stars, cannot be a free gas 
of quarks (jOzel et al.l l2010l : IWeissenborn et al.l l201ll ) , an intera cting phase of quarks still appears 
to be consistent with the recent finding of a 2Mq neutron star ( Demorest et al.ll2010l ). 



The main questions we seek to answer are: given a few initial parameters of the newly born 
neutron star, such as its spin period, magnetic field and baryonic mass, can we determine which 
neutron stars are likely to eventually manifest a quark phase in their interior? If so, how long 
before the transition to quark matter occurs? In this work, we take a step towards answering 
these questions by taking a closer look at the rotational evolution of a newly-born hot neutron star 
as it spins down, and we focus on two main driving factors - magnetic braking and gravitational 
radiation from r-modes. We consider different equations of state (EoS) for neutron stars, since 
the EoS at high density is uncertain and lacking strong empirical constraints (although a recent 
statistical analysis shows that the equation of state s tiffens at high den sities and is consistent 
with the expected range in certain nuclear parameters: ISteiner et al.ll2010l ). Furthermore, second 
generation gravitational wave detectors such as Advance d LIGQ will soon be operational, and our 
work provides an update for a similar theoretical study (JHo &: Laill2000l ) performed several years 
ago with regard to the LIGO detector. We should note here that in contrast to the assumptions 



in IHo fc Lail ([20001) , our work includes the effects of magnetic damping of r- modes in the spindown 
evolution, which leads to quantitatively different results. 

In section [21 we outline the approach and main equations tha .t describe theneutron star's 
spindown. The theoret ical analysis follows in pa rt the works of [Ho &: Lai[ (|2000l ) , as well as 



Cuofano fc Dragol (|20ld ) and [Rezzolla et al.[ (|2001bl ). Section [3[ collects our main conclusions from 
this analysis. In section [H we analyze the e volution of the gra vitational wave frequency associ- 
ated to the growing r-mode, along the lines of [Owen et al.l (|l998l ). and present results in section [5l 
We conclude in section [6[ with a preliminary calculation of the gravitational wave signal from an 
explosive deconfinement transition (the "Quark-Nova"). 



2. Spindown to deconfinement 



In a previous paper (IStaff et al.[ [20061 ). we had determined birth parameters of a neutron star 
which would support the transition to deconfinement driven solely by magnetic braking. We found 
that neutron stars that are born with mass M > 1.5Mq and spin-period -P ^ 3 ms are the best 
candidates to reach deconfinement density in their core (assumed to be ~ S/jq), given a range 
of neutron star magnetic fields 10^^-10^^ Gauss. This small value of the initial spin-period (or 
large spin frequency) is required since magnetic braking from a rotating magnetic dipole is itself 
proportional to the third power of the spin frequency. Most of the increase in central density occurs 
soon after the rapidly rotating neutron star is born. 

In that work, we neglected the role of gravitational radiation from r-modes in an adhoc way - 
by assuming that neutrons stars are born axisymmetric and remain so. In reality, non-axisymmetric 
perturbations are expected due to the violeiit pro c ess of a supernova, which leaves the new-born 
neutron star in a turbulent state (|Keil et al.l [19961 : [Akivama &: Wheeled [20061 ) . What is then the 
role of r-modes and gravitational braking in comparison to magnetic braking as far as the time 
to deconfinement is concerned? Previous studies on the interplay of r-modes and neutron star 
magnetic fields have focused on different questions - for example, [Led (l2009l ) has studied the r- 
mode in magnetized neutron stars with B < 10^ ^ G with a r i airn to understand X-ray pulsations 
from local hot spots in accretion-powered pulsars. [Ho &: Lail (120001 ) have shown that magnetic fields 
oi B > lO^^G can make magnetic braking as important as r-mode spindown, specially for slowly 
rotating stars, and that Alfven wave driving of the r-mode can also play a role. 



Meanwhile, other works ( Rezzolla et al.l boOlal lbl: iKiuchi et al.ll201ll : ICuofano &: Dragol I2OI0I ) 
have focused on the evolution of toroidal magnetic fields generated by the secular r-mode. Essen- 
tially, the toroidal field is generated by differential rotation of the stellar fiuid associated with the 
r-mode. The differential rot at ion is a key fea t ure of the r-mode instability in the non-linear regime 



( Levin &: Ushomirsky 



2001 



Rezzolla et al.l ([2001al ) have shown that, for isolated neutron stars. 



this effect can amplify existing magnetic fields by two orders of magnitude within the time taken 
by the r-mode to saturate (w few hundred sees.). Back-reaction on the r-mode due to this toroidal 
field implies that there is an associated magnetic damping, which has to be factored into the spin 



evolution. For neutron stars that are accreting from a b inary companion, the bu ild-up of toroidal 
fields to sizeable values is rather slow, and according to ICuofano &: Dragd (J20ld ). can take up to 
several hundred years. However, we are interested in this effect on isolated neutron stars. The 
magni tude of the magnetic damping term Fm is proportional to the integrated time evolution of 
a'^{t) ( Rezzolla et al.ll2001bl ). with a being the r-mode amplitude. These authors obtained a strong 
damping effect by assuming that the background evolution of a{t) is identical to that given by the 
absence of any magnetic fields. However, we solve the evolution equations for the star's angular 
frequency fi and a{t) with magnetic fields and magnetic damping from the start. Fm grows while 
the mode is unstable, and once it saturates Fm remains constant. This consistent inclusion of F^ 
leads to a slower evolutionary path for a{t) towards saturation, hence the magnetic damping is not 
as effective in suppressing the r-mode. Consequently, for a range of magnetic fields and small r- 
mode saturation values, we find that this timescale can be much larger than the typical time taken 
by the neutron star to spindown to typical quark deconfinement density. Therefore, our results 
clearly indicate that the r-mode, for realistic saturation values, quantitatively affects the spindown 
evolution of an isolated neutron star. We now discuss the relevant equations for the evolution of 
the star's rotation rate and r-mode evolution. 

The spindown fi of a neutron star (mass M, radius R) is accompanied by a loss of energy 
as well as angular momentum. For the radiating star, conservation of total angular momentum 
Jtot = -^tar + {^ — Kj)Jc with Kj a dime nsionless constant and Jsi,a,r = I^ (where / is the moment 



of inertia about the rotation axis) yields (IWagoner 



2002 



Cuofano fc Dragoll2010l ) 



dJ, 



tot 



dt 



2JcF„ + Ja- I^F, 



mag ■ 



(1) 



Jc = —Kco'^Jsta.r is the canonical angular momentum of the r-mode to 1st order in Q, (JFriedman fc Schutz 



19781 ). Fg is the rate of gravitational radia tion associated to the l =m=2 current multipole taken 



for an n=l polytropic star from eqn.(65) of lAndersson et al.l (|200ll ). This is expected to be a good 
approximation to a large part of the neutron star interior. J a is the accretion rate onto the star 
(assu med zero for our case since we consider isolated stars only) and i^mag is the magnetic braking 
rate (JManchester &: Tavlorill977l ). The dimensionless quantity K^ = 3 J/2/ is defined from 



J 



1 



R 



MR"^ Jo 



pr dr . 



pr dr . 



(2) 
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where p = p{r) is the density profile (assumed radially symmetric) of the star Q|. J^ evolves 
in time as a result of competing influences from gravitational damping (which feeds the r-mode), 



^ In principle, equatorial flattening due to rotation renders the density proflle asymmetric, but this does not change 
the value of Kc significantly, which is the basis for our approximation of a symmetric profile. 



viscosity (which damps the r-mode) and the magnetic damping term Fm- Both bulk and shear 
viscosities are included in our analysis although for simplicity we keep the star at a uniform tem- 
perature of lO^K for the duration of the evolution. Although viscosities (especially bulk viscosity) 
is strongly temperature-dependent, our approximation is not as drastic as it may seem. A typical 
cooling p rofile of a neutron star, driven by neutrino emission from the modified URCA process is 
given by (|Qwen et al.lll998l ) 
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+ T, 
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(3) 



where Tg is the star's core temperature T in units of lO^K, Tj^g the initial temperature, and 
Tc ~ 1 yr is the characteristic cooling time from the modified URCA process. If we begin with birth 
temperatures T ~ lO^^K, we see that within a few seconds, we have T ~ lO^K. Since the r-mode 
does not have a large impact until t > 10^ seconds, we can approximate Tg = T/(10^K) ~ 1. 
Subsequent cooling is on the timescale of years a nd is also a small pertu rbation on our results. 
Using the evolution equation for Jg, (see eqn.(4) of lCuofano k. Dragd 120101 ) . one can write 
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2a'' K, [KjFg + (1 - K,)[F, + F„ 
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Fmag and 
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a ^ [^^-[^^ + ^-]]-^' 



where F„ is the viscous damping rate (JAndersson et al.M200lh. We solve the two coupled 
equations above numerically. We use the RNS code (jStergioulas &: Friedmanlll995l ) to construct 
2-dimensional models of rapidly rotating neutron stars. For a given EoS and for a fixed baryonic 
mass, the RNS code outputs a sequence of neutron star models (with specific gravitational mass, 
radius, spin frequency etc.) that have increasing central density and decreasing angular velocity. 
The fastest model in such a sequence spins at or near Kepler frequency. Note that the magnetic 
field does not appear in the RNS code, and in any case, its effect on structure at the field values 
considered here is negligible. The magnetic field only de termines th e time to deconfinement, unless 
the spindown is completely r-mode dominated. As in IStaff et al.l (J2006l ) we assume 5po to be a 
critical density at which quarks in the interior of the neutron star become deconfined - this is not a 
well determined n umber, and rnay sp an a range from (4-8)/9o if a mixed phase of quark and nuclear 
matter is favored (jGlendennine] 1 19921 ) . To counter this uncertainty, we have checked our numerical 
results for a higher putative deconfinement density {p ~ 8po) and found that it does not change 
any of our quantitative conclusions by a significant amount. 

Using the RNS code, we construct sequences of stars (for a given EoS) with constant baryonic 
mass and decreasing spin such that the non-rotating configuration has a central density equal to 
5po- We then calculate the time the rotating star takes until its central density is within 99% of 
5/90- This is a practical way of using eqn.(j3|) to obtain the minimum mass required (for a given 
EoS) to support deconfinement, since the time taken to reach zero spin is infinite from the magnetic 



braking being proportional to 0. Since we are interested in the time it takes for a star to spindown 
to the critical density for deconfinement, by choosing a sequence where the non-rotating model has 
a central density equal to the critical density, this time gives us the maximum time the star takes 
to reach the deconfinement density. 

We assume a very small initial amplitude of the r-mode for the fastest spinning star in a 
sequence, with typical initial a ~ 10~^, although our results for those cases where the r-mode 
saturates are insensitive to this initial value. For the case when the mode does not saturate, we 
do find mild sensitivity to the initial value of a. From the first of eqn.Q, we can calculate the 
time-step between two consecutive neutron star models in the spin-sequence output by RNS code 
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-2K,n{U)a{UY[K,Fg + (1 - Kj){F, + Fm)] - n{U)F„ 



(5) 



The inverse dependence on a implies that for small magnetic fields, where -Fmag is small, At is 
large when the mode amplitude is just starting to grow. The RNS code then needs to generate two 
widely separated rotating configurations within a sequence, leading to low resolution in some parts 
of the P vs t curves for small magnetic fields (eg., the curve for B=W^^ G in FigH]). However, 
this is only a problem initially and does not affect the total time to deconfinement. The above 
equation, along with the second of eqn.(jl|) is then used to determine the time evolution of a. For 
our analysis we have used four different magnetic field strengths equally spaced on a logarithmic 
scale from 10^^ G to 10^^ G. We display results for three different equations of state: EoS A which 
is composed of only neutrons (|Pandharipanddll97lh and uses a variational principle to determine 



the minim um energy state; EoS BBB2 which includes muons and uses a field-theoretic many-body 
approach (JBaldo et al.lll997l ) ; and finally EoS APR wh ich is relatively stiff and admits a mixed phase 
of quarks and nuclear matter for the heaviest stars (lAkmal et al.lll998l ). All three EoS generate 
stable configurations for central densities exceeding 5po and can go up to even 8po. They differ in 
the details of the density profile, with the softest (EoS A) having higher interior densities for the 
same mass (more c ompact). The maximum gravitational mass exceeds 2 M(d only for th e APR EoS 
(JAkmal et al.lll998l ). with BBB2 providing a maximum mass of 1.92M0 (jHaenseJ I2OO3I ) . EoS A is 
too soft to generate a 2M0 n eutron star and may s eem inadequate to explain recent observations of 
such massive compact stars (JDemorest et al.ll201Cll ). However, if such massive stars are really quark 
or hybrid stars, we cannot rule out EoS A in this way, and its inclusion is still useful simply to 
examine the trend of a relatively soft EoS on the gravitational signal from r-mode driven spindown. 



3. Results: r-mode spindown versus magnetic braking 

The following results are obtained on solving eqns.(4) numerically. In Fig. [H the curves 
show how the period, starting from the Kepler rotation rate at birth, evolves as a function of 
magnetic field strength, for different equations of state. Note that these curves correspond to 
structural parameters of that particular stellar configuration which just reaches the putative quark 
deconfinement density at zero angular velocity - we may call this the critical configuration. For 



6 



any given EoS, this is how we determine the minimum n iass and spin peri od required for a quark 
phase to appear as a result of spindown. UnHke Fig. 4 of IStaff et al.1 (J2006l ). spindown is no longer 
always magnetic-dominated for arbitrary magnetic fields. For example, for B < 10^^ G, r-mode 
spindown takes over at a few hundred seconds. As the instability develops, magnetic damping Fm 
grows and limits the growth of the r-mode. This leads to a plateau in the period that can last up to 
several years before magnetic damping once again becomes the dominant driver of spindown. It is 
noteworthy that for the critical configuration, there is no change in the time to deconfinement for 
any value of B. However, if we choose a heavier mass at birth, the time to deconfinement would be 
shorter, and can occur while the star is still being spun down as a result of the r-mode. Assuming 
that even the most massive neutron stars can be transformed to stable quark stars if they reach the 
deconfinement density (that is, the quark matter EoS should be sufficiently stiff to support such a 
mass), we may expect that several such neutron stars would have already spun down to the point 
where they underwent a quark-hadron phase transition, and now contain quark matter in their 
core. 



Our findings here are different from the conclusions of lHo &: Lail (|200Cll ) who studied the r-mode 
of magnetized and slowly rotating neutron stars, and found that the r-mode alters the spindown 
if magnetic fields are less than B ~ lO^^G. Due to the inclusion of the magnetic damping term, 
which was omitted in their work, we find that r-mode driven effects on spindown are pronounced 
only for B < lO^^G, since for higher B fields, magnetic damping effectively kills the r-mode. We 
emphasize that the effects of magnetic da i nping f or smaller B fields are not as severe as might be 
expected from the work of iRezzolla et al.l ()2001al ) due to the self-consistent evolution of a in the 
presence of the magnetic damping term, as explained at the beginning of section [2j 

We also find that the r-mode driven spindown is more pronounced in a relatively stiff equation 
of state such as EoS APR. This is because Fg, the gravitational damping rate is larger for a stiffer 
EoS G, driving the r-mode unstable in a shorter amount of time (effectively spinning the star down 
to deconfinement density quickly). 

The saturation amplitude asat of the r-mode assumed in the foregoing analysis is 0(1). Such 
a large value is p r obably unrealistic due to multi-mode coupling and onset of non-linear effects 
(jArras et al.l 120021 : iBondarescu et al.ll2008l ). If we use a smaller saturation value, we find that our 
conclusions change quantitatively. The r-mode grows more slowly for a smaller asat; but continues 
to be important since the main damping agent Fm, which is roughly proportional to a^, is also 
smaller. As a consequence, the spindown due to the r-mode is weaker, but also lasts longer, as 
seen in FigiS) This can lead to the star taking a shorter time to reach a given period in the case of 
smaller B fields. This can be seen clearly in Figl6] which compares the time taken by the star to 
spindown to a period P=3 ms, for asat=0.01 and asat=0.005. 



^ Fg oc MR^ /P^ and for stars with the same baryonic mass and central density at zero angular frequency, this 
quantity is systematically larger for a stiffer EoS. In other words, for a given baryonic mass, a stiffer EoS can support 
a larger gravitational mass and radius, and has a slightly smaller Kepler frequency since its average density is smaller 



APR / niEoSA 



(less compact). For eg., Fg /F, 
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4. Evolution of the r-mode and Gravitational Wave Frequency 



As shown in the previous section, fairly large magnetic fields of 10^^ G or more are required 
to make the r-mode irrelevant for spindown. For smaller magnetic field s, the r-mod e evolu tion can 
be the main agent for spinning down rapidly rotating stars. Following lOwen et alj (|l998l ). we can 
obtain expressions for the growth of the mode amplitude (a) and the corresponding evolution of 
the frequency of the gravitational wave (/) under certain approximations. 



4.1. Saturated Phase 



The r-mode ampli tude grows rapidly until non-linear saturation occurs (jArras et al.l l2002l : 
Bondarescu et al.l 120071 ) at which point d=0. We can obtain a description of this saturated phase 
from eqns.dH), which imply that the rotation rate evolves according to 



n 



2a'K,F„ 



n l-a^Kc{l-Kj)' 



(6) 



where a = Ogat is the saturation amplitude of the r-mode. Note that Fm does not appear since it 
can be eliminated using the second line of eqns.(4). Its effect will however show up in the growth 
phase discussed below. Using the fact that Q, = 37r//2 for the I = m = 2 r-mode, setting Kc = O.l 
(a typical value determined from eqn.(l2])), and using the frequency-dependent expression for the 
gravitational wave damping timescale tgr = -^g~^! '^^ find 
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where we have used the expression for tgr from lLindblom et al.l ( 19981 ) and neglected a^Kdl — 
Kj) <^ 1 in the denominator of eqn.Q. B, R and I are expressed in terms of reduced dimensionless 
units. 



4.2. Growth Phase 

In the initial stages of the mode growth phase, the viscosity controls the evolution of the 
r-mode, but the instability to gravitational waves soon takes over. Viscous damping from bulk 
viscosity can be large at T > lO^^K, but as the neutron star cools rapidly on the order of seconds, 
our approximation of setting Tg = 1 implies that viscosity does not affect the subsequent spindown 
behaviour or the gravitational waveform during the growth of the r-mode instability. It follows 
from eqns.(j4]) that during this phase, since a is small 



^-(F -F )- — 



(8) 



The angular velocity evolves according to the first line of eqns.(4). We assume Kj « for 
the growth phase, which is tantamount to including fully the canonical angular momentum of the 
r-mod e in the star's ph ysical angular momentum. This is only justified if differential rotation is 
small (ISa &: Tomdl2005l ). as is the case when the mode is still small but growing. Then, as before, 
the proportionality between / and U implies 
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where we have used the expression for Fm from eqn.(ll) of ICuofano &: Dragd (|20ld ). We 
can now use eqns.([7I) and ([9]) to obtain the gravitational strain amplitudes and waveforms for the 
saturated phase and growth phase respectively. 



5. Results: Gravitational Strain Amplitudes and Waveforms 

The strain amplitu de h(t) corresponding to the l=m=2 mode is found by a standard multipole 
analysis (|Thorndll98d ). Including the angle-average over the position of sources in the sky, 
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(10) 
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where ^22 is the current multipole as given by eqn.(3.9) of lOwen et al.l (Il998l ). the mode 
frequency w = 40/3 and D is the source distance, chosen henceforth to have a typical value of 20 
Mpc|j. Figure [2] shows the time-evolution of the strain amplitude h{t) for various magnetic fields 
and for different EoS for Osat = 0.5, while Fig. [7| shows this for agat = 0.01. For magnetic fields 
B > lO^^G, the r-mode withers much before it nears saturation, for any EoS - the strain amplitude 
is consequently very small. For B ~ 10^^ — lO^^G, we notice a strong dependence on the equation 
of state. The soft EoS A does not lead to a saturating r-mode while relatively stifFer EoS do. For 
weaker mag netic fields around B 
in Fig. 5 of lOwen et aP (|l998l ). 



10 G, the r-mode saturates and displays the behaviour shown 



We now examine the gravitational waveform in the frequency domain h{f) which is the Fourier 
transform of the time signal. This is useful in estimating the signal-to-noise ratio (SNR) in matched 
filtering techniques. In the stationary phase approximation. 



^This fiducial value of D extending t o the Virgo cluster is chosen to encompass enough neutron stars to ensure a 
reasonable event rate I Owen et al.lll998l '). 
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From eqns.d?]), ([9]) and (jlOp . we can find h{f) for botli ttie growtfi pliase and the saturated 
phase. From Figs. [3] and [8] (for Ogat = 0.5 and 0.01 respectively), we see that for magnetic fields 
B < 10^'^G, th e signal frorn. the r- mode is similar to previous result obtained in the absence of 
magnetic fields (lOwen et al.lll998l ). A sharp peak at high frequency (~ 1.6 kHz) corresponding to 
the growth phase is seen, followed by a plateau signal at lower frequencies as the r-mode saturates. 
Note that the curves terminate on the left since at that point, deconfinement density is reached 
and the subsequent signal from the phase transition needs a detailed analysis beyond the scope 
of this work. Howeve r, we refer the reader to the Appendix for a preliminary treatment of this 



issue. Previous works (jSchenk et al.ll2002l : iMorsinkI 



2OO2I : 



Brink et al.ll2004l : iBondarescu et al.ll2007l ) 



indicate that the saturation amplitude Ogat may be quite small, around 10 ^ instead of order 1. 
Still, it is useful to contrast both cases: when Osat ~ 0.5 and asat ~ 0.01. 



5.1. Signal Detectability 

In Fig. m we present our estimate of the continuous gravitational wave signal from the spin- 
down (assuming a point source at a distance of 20 Mpc), compar ed to the anticipated sensitivity 
of Advanced LIGO|j and the Einstein telescope ( Hild et al.l 120101 ). The left panel of Fig. [5] shows 
the weighted strain amplitude (1/vHz) versus frequency for the three EoS studied in this paper 
at 10^^ G, compared to the anticipated noise- weighted design sensitivities of Advanced LIGO and 
the Einstein telescope. The right panel shows the same for 10^"^ G. In both figures, we are assum- 
ing Qfsat ~ 0.01. The saturated phase of the curve is at least an order of magnitude above the 
anticipated sensitivity for Advanced LIGO for all EoS at 10^^ G and could therefore be expected 
to be observed out to 20 Mpc with Advanced LIGO, with a typical SNR of ~ 20. At 10^^ G and 
for the softest EoS, the SNR in Advanced LIGO drops to ~ 2. The Einstein telescope is expected 
to gain about an order of magnitude in sensitivity compared to Advanced LIGO, and hence all of 
the above cases should be detectable in principle. Our SNR estimates only provide upper limits 
since matched filtering is probably impractical for such sources, where spin parameters cannot be 
promptly measured, though the source location could be known if the supernova is observed. How- 
ever, we note that most of the signal comes from the saturated phase, which for the case of small 
asat occurs several years after a neutron star's birth in a supernova. Observations of the neutron 
star's spindown parameters may then be possible, making the matched filtering method more fea- 
sible. In the context of r-modes in newly-born neutron stars, methods other than matched filtering 



^These curves represent the incoherent sum of the principal noise sources, such as quantum noise, seismic noise and 
thermal noise, as best understood at this time. There will be, in addition, technical noise sources. These curves serve 



as a g uide to the overall curve and an early approximation to the anticipated sensitivity (JLIGO Scientific Collaboration 

horn . 
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have been suggested (JBradv &: Creightonlll998l : lOwen et al.1 l20ld : IZhu et al.ll201ll ) that would be 
less computationally intensive. 



5.2. Saturation at Osat ~ 0.5 



As mentioned in I Owen et al.l (|l998l ) , the sharp spike in the strain amplitude seen in certain 
cases only lasts for a brief period of time (of the order of a minute) and carries a small fraction of the 
total energy emitted in gravitational waves over the spindown epoch (see Fig. [2]). However, for the 
softer EoS, such as BBB2, we see that for magnetic fields B ~ lO^^G, the spike is softened into an 
asymmetric hump. In this case, the growth phase lasts about an hour and carries a larger fraction 
of the emitted energy - it may be detectable owing to its larger SNR in second and third generation 
detectors. Furthermore, for the softest EoS considered here (EoS A), there is no saturation regime 
for fields larger than B ~ 10^^ G. These are potentially distinguishing feature of the equation of 
state in gravitational waves. In general, for very high magnetic fields B > 10^^ G, the r-mode is 
strongly suppressed, and the signal in gravitational waves is weak, with little possibility of detection 
even in third generation detectors. 



5.3. Saturation at asat ~ 0.01 

For a more realistic value of Osati we find that the strain amplitude is about an order of magni- 
tude smaller than for asat=0.5, but the signal persists for a much longer time (several months/years 
for B < 10^^ G). This is clear from Fig. [71 This is because of the dependence of Fm on B^ and 
a{t). For smaller B fields and asat, the weaker magnetic damping allows the r-mode to survive for 
a longer time, even though mode growth is slower. This eventually leads to a faster spindown for 
the star, compared to a higher Ogat- To emphasize this effect, in Fig. [U we compare spindown with 
asat=0.01 and 0.005, for the same initial stellar configuration and the same magnetic field. The 
crossing of the two curves in that figure illustrates this effect. However, we note from comparing 
Figs. [3] and [8] that h remains relatively independent of asat, although both the B = 10^'^ G and 
B = 10^^ G cases shows a flat segment corresponding to the saturated phase for Ogat = 0.01. Es- 
sentially, this is because the decrease in the peak value of h{t) for smaller asat is compensated by 
the corresponding decrease in / (eg., as in eqn.([7])). This implies that we require a larger interval of 
time integration to validate the stationary phase approximation for a smaller asat (eqn. [T2[). Thus, 
h remains almost unchanged. 



6. Conclusions 

We have studied the role of r-modes in spinning a rapidly rotating and magnetized neutron 
star down to typical quark deconfinement densities. Apart from the usual spindown associated to 
magnetic braking, the magnetic damping of r-modes plays an important role in the period evolution 
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of the star. We find that, for reahstic (small) values of the r-mode saturation amplitude Osatj the 
time to deconfinement is sped up by the r-mode instability when magnetic fields are of order 10^^ G 
or less. This result reflects the strong dependenc e of the inagnet ic damping effect on B and the 
evolution of a. Thus, in contrast to the result in (JHo &: Laill2000l ). where the r-mode was seen to 
affect the spindown already for magnetic fields less than B ~ lO^^G, the inclusion of the magnetic 
damping term leads us to conclude that r-mode driven effects on spindown are pronounced only 
for B < lO^^G. We also explored the gravitational wave signal generated during the spindown 
phase as the r-mode first grows then saturates - we follow the signal continuously until the quark 
deconfinement threshold is reached. For realistic values of agat, the r-mode saturates and leads 
to a strong signal in gravitational waves, except in case of very soft equations of state and very 
large magnetic field (10^^ G or more). There could be a sizeable fraction of the total rotational 
kinetic energy emitted as gravitational waves in this epoch, which can last several years, making 
it detectable in upcoming second and third generation detectors. Therefore, gravitational waves 
could be used in this manner to probe the equation of state inside neutron stars. As shown here, 
using r-modes is an alternate way in which the EoS can be probed through gravitational waves u 

Gravitational waves are going to open a new window of observation into our universe. Among 
the many discoveries that will be made, we anticipate the exciting prospect that neutron star 
spindown will reveal signs of the elusive r-mode instability as well as signatures for the onset of 
quark matter in the core. In the appendix we outline a first preliminary estimate of what the 
gravitational wave signal from a "Quark-Nova" would look like. It would be interesting to examine 
this signal with detailed simulations, especially in cases when the Quark-Nova occurs very shortly 
after the neutron star is born in a Supernova - the signature of this " dual-explosion" in gravitational 
wave detectors would be two very different signals coming from the same source but separated in 
time by a few days to few weeks, depending on the time delay between the two explosions. Such a 
signature would be unmistakeable in upcoming gravitational wave detectors. 

Appendix: Gravitational waves from the quark-hadron phase transition 

It remains an open question as to what happens to the neutron star when the deconfinement 
density is reached. One possibility is that the entire sta r is converted to a quark star due to the 
inherent stability of strange quark ma tter (IWittenI lOS^) . If this conversion occurs in an explosive 
manner, a "Quark-Nova" could result (jOuved et al.ll2002l ) - what will the gravitational signal from 
such an event look like? The conversion involves a two-stage process - neutrons in the core dissolve 
into a u and d quark fluid that is more compact, causing the core to shrink, followed b y combustion 
to u, d, s quarks through leptonic as well as non-leptonic processes. iLin et al.l (J2006l ) have studied 



^We note that recently derived empirical ugisej: bounds on the R ravitation al pow er radiated by the Crab pulsar 
constrain the ellipticity of deformed pulsars (iPitkinil2011 : Abbott et al. 2008j, 20101 ) with implications for the max- 
imum theoretical elasticity of the neutron star crust ( Horowitz fc Kadaul 120091 ). On the other hand, in order to 
similarly constrain the r-mode amplit ude frorn sear ches directed at known pulsars, a different range of frequencies 
and polarizations must be probed (jOwen et al.ll2010l ). thus the limits placed from the Crab pulsar do not have direct 
bearing on r-modes. 
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the first stage with Newtonian hydrodynamics and found that quadrupolar and quasi-radial modes 
are excited during the collapse, leading to gravitational wave emission with an energy output of 
~ 10^^ ergs. However, this work applies to a mixed phase of quarks and nuclear matter and does 
not consider an explosive transition that begins with non-premixed fluids. 



A full numerical treatment of the second stage, an explosive phase transition taking into ac- 
count fluid motion in 3D is a complex tas k and beyond th e scope of this work. However, preliminary 
steps have been taken in this direction. iNiebergal et al.l (|20ld ) solved hydrodynamical flow equa- 
tions for the combustion of neutron matter to strange quark matter in the laminar approximation, 
including weak equilibrating reactions and strange quark diffusion across the burning front. The 
numerical results suggest laminar speeds of 0.002 — 0.04 times the speed of light, muc h faster than 
previous estimates derived using only a reactive-di ffusive descriptioii (IQliiitd 119871 ) . Turbulent 
combustion has been addressed in a recent work bv lHerzog &: Roepkd (|201ll ) who found that the 
combustion stops short of converting the e ntire star to quark m atter (the reaction is no longer 
exothermic) . This was also the conclusion in INiebergal et al.l (120101 ) , though for a different physical 
reason |2|. Neither of these works continue on to estimate the gravitational signal from the explo- 
sive combustion. Following the hydrodynamical approach of these studies, we have estimated the 
gravitational wave signal h{t) from the following steps. 



Starting from a mechanically stable configuration where quark matter con stitutes a small 



fraction of the stellar core, we initiate combustion at the speeds obtained in INiebergal et al. 



(|2010l ). The quark fluid, described by a simple bag model equation of state is subsequently 
evolved using inviscid hydrodynamical equations for relativistic fluid flow (relativistic Euler 
equations) coupled with Newtonian gravity |I|, assuming an axisymmetric rotating configura- 
tion of the star about the i-axis. In the cylindrical coordinates (r, (/>, 0), where (p is the polar 
angle and 9 is the azimuthal angle, the relevant equations are given by dt{rU) + dr{rF) + 
dz{rG) = S , where 
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and where D = pj, S = D/17, and r = S—T)—p are introduced solely for the purpose of writing 
the relativistic Euler equations in a form analagous to the more familiar non-relativistic Euler 



iNiebergal et al.l ( 201(1 ) found that as the burning front expands and cools, it enters an advection dominated 
regime, where the upstream (hadronic) fluid velocity advects the interface backwards faster than it can progress due 
to reactions and diffusion. Consequently, the interface halts short of the neutron star surface. 

^ Although the Poisson equation for gravity violates the speed of light (since it is an elliptic PDE), it is unavoidable 
unless one uses the full GR. 
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1, the fluid velocity, Lorentz factor, gravitational 

{Ur,Ue,Uz), 7 



equations. In the natural units where c 

vector field, and specific enthalpy are respectively given hy u = {ur, ug, u^), 7 = l/y^l — jup, 
g = — V</>, and h = 1 + e + p/p, where e is the specific internal energy, p is the pressure, 
and (p is the gravitational potential. In the non-relativistic limit (7 — )• 1), we recover the 
non-relativistic Euler equations. 

The quark matter has a density pq = bpn and is described by the Bag model (bag constant 
B^i'^ = 145 Mev) while the nuclear fluid has density pn = Ps&t = 2.5 x lO^^g/cm^ and is 
described by the perfect caloric EoS with an adiabatic index of 1.7. At t=0, we initiate 
combustion with a density discontinuity pq — pn, which is situated at a radial coordinate 
r = R/A inside the star, and choose an initial burning front sp eed of Ur-=0.01c, which is in 
the range of burn velocities obtained in iNiebergal et al.l (120101 ). We then evolve the above 
equations using a weighted average flux relativistic HLL (Harten, Lax, and Van Leer) solver 
on a cartesian grid with spatial resolution 0.5 km and a timestep determined by a Courant 
number of 0.3. The solver is coupled with the Poisson equation V'^cf) = AirGp. 



Solving for u{r, z, t) and p{r, z, t), we follow I Zwerger fc Mullerl (jl997l ) to compute the quadrupole 
wave amplitude Aon( t) for our axisymraetric configuration, and from there the signal h{t) as 
given by eqn.(21) of (IZwerger &: Mulleiill997l ). 



^ ^ SV vr D 



(1) 



We assume a source distance D = 20 Mpc to plot the maximal signal strength h(t) and the 
corresponding luminosity L{t) in Figl9l The origin of the peak in h{t) is a balance between 
increasing mass outflow and a decreasing density gradient. The luminosity is proportional to 
the square of h{t) (hence the sharp dip at the maximum of h{t)) while the total emitted energy 
obtained from integrating the luminosity curve is ~ 1.4 x 10^* ergs, which is about 5 orders 
of magnitude smaller than the binding energy of the neutron star (10^'^ ergs) and 2 orders of 
magnitude less than the energy released in gravitational waves during the first sta ge of core- 



collapse, where the energy coin es from coupling oscillations to rotational motion (JLin et al. 
200d : iMarranghello et al]l2002l ). 



Based on this result, the gravitational wave signal for this stage of the phase conversion would 
be hard to detect in Advanced LIGO, unless the source is Galactic (within few kpc). Most of the 
energy released in the phase transition is in the form of latent heat and neutrinos. However, our 
results are only a preliminary estimate, designed to provide a guide to the expected signal from a 
Quark-Nova. Simulations of gravitational wave signals from realistic core-collap se supernova models 
indicate that convective flows driven by neutrino heating (JMiiller et al.ll2004l ) can drive a strong 
gravitational wave si gnal. The typical tim e for neutrinos to diffuse out of the hot quark star (~ 
0.1-1 sec according to lKeranen et al.ll2004l ) is about an order of magnitude larger than the duration 
of the signal due to explosive combustion to strange quark matter, therefore, we can expect neutrino 
heating to be important in our context as well. We have not included such convective effects in 
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our present simulation, nor the effect of magnetic fields, which can also modify Rayleigh- Taylor 



instabilities (JLugones et al.ll2002l ). 
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Fig. 1.— Period vs time for the APR EoS {top left), EoS A {top right), and EoS BBB2 {bottom left) 
with a constant temperature of T = 10^ K. The curves from left to right are labelled for log(i3). 
For B > lO^^G, magnetic braking starts to dominate the r-mode as far as spindown is concerned. 
The curves terminate abruptly at the moment when the central density is within 1% of the critical 
density (see discussion above eqn.([5])). 
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Fig. 2. — Gravitational wave strain as a function of time for the APR EoS {top left), EoS A (top 
right), and EoS BBB2 [bottom left) with a constant temperature of T = lO^K. The curves for softer 
EoS and larger B fields have a much lower peak strain as in these cases the r-mode amplitude never 
saturates. 
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Fig. 3. — Gravitational wave strain as a function of frequency for the APR EoS (top left), EoS A 
{top right), and EoS BBB2 {bottom left) with a constant temperature of T = lO^K. The 10^^ G and 
10^^ G curves are distinct from the rest, as in these cases the r-mode amphtude barely or never 
saturates. 
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Fig. 4. — Gravitational wave strain as a function of frequency for three different EoS - APR(solid), 
BBB2 (dash-dotted) and EoS A(dotted), compa red against sensitivity curves (strain noise) for (grey) 
Advanced LIGO taken from Harrv et a.1. (120101 ^ for the LIGO Scientific Collaboration, and (black) 
the Einstein Telescope from lHild et al.l ( 2010l ). Left panel is with B=10^^ G and right panel with 
B=10^^ G. The saturation amplitude Q:sat=0.01. 
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Fig. 5. — Same as Fig. [H but with r-mode saturation value Osat = 0.01. Note that the effect of 
the r-mode is still significant but appears later than for asat = 0.5. This is due to the nature of 
magnetic damping. 
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Fig. 6. — Spindown from Kepler frequency for identical stellar configurations (APR EoS with 
B=W^'^ G) but with different r-mode saturation values: asat=0.01 and 0.005. Time to reach P=3 
ms (where supposed deconfinement density 5po is reached) is shorter for smaller asat- 
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Fig. 7. — Same as Fig. [2] but with Osat = 0.01. For smaller values of 5, the signal reflects the fact 
that the r-mode evolves on much longer timescales than for Ogat = 0.5. 
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Fig. 8.— Same as Fig. El but with a^at = 0.01. The 10^^ G (and for EoS A also the 10^^ G) are 
distinguished, as in these cases the r-mode amplitude never saturates even when we assume such a 
low saturation value. 
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Fig. 9. — Gravitational wave strain h{t) and luminosity L(t) as a function of time for combustion 
of neutron matter to strange quark matter following quark deconfinement inside a neutron star. 
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